Enhancing prognostic power in multiple myeloma using a plasma cell signature derived from single-cell RNA sequencing.

Multiple myeloma (MM) is a heterogenous plasma cell malignancy, for which the established prognostic models exhibit limitations in capturing the full spectrum of outcome variability. Leveraging single-cell RNA-sequencing data, we developed a novel plasma cell gene signature. We evaluated and validated the associations of the resulting plasma cell malignancy (PBM) score with disease state, progression and clinical outcomes using data from five independent myeloma studies consisting of 2115 samples (1978 MM, 65 monoclonal gammopathy of undetermined significance, 35 smoldering MM, and 37 healthy controls). Overall, a higher PBM score was significantly associated with a more advanced stage within the spectrum of plasma cell dyscrasias (all p < 0.05) and a shorter overall survival in MM (hazard ratio, HR = 1.72; p < 0.001). Notably, the prognostic effect of the PBM score was independent of the International Staging System (ISS) and Revised ISS (R-ISS). The downstream analysis further linked higher PBM scores with the presence of cytogenetic abnormalities, TP53 mutations, and compositional changes in the myeloma tumor immune microenvironment. Our integrated analyses suggest the PBM score may provide an opportunity for refining risk stratification and guide decisions on therapeutic approaches to MM.

Donor germ-line variants associate with outcomes of allogeneic hematopoietic stem cell transplantation in patients with myelodysplastic syndromes.

Overall survival probability for MDS patients who underwent allo-HCT and were matched to donors that are wild-type (red) and heterozygous (blue) for the rs111224634 SNP.

Improving polygenic risk prediction in admixed populations by explicitly modeling ancestral-differential effects via GAUDI.

Polygenic risk scores (PRS) have shown successes in clinics, but most PRS methods focus only on participants with distinct primary continental ancestry without accommodating recently-admixed individuals with mosaic continental ancestry backgrounds for different segments of their genomes. Here, we develop GAUDI, a novel penalized-regression-based method specifically designed for admixed individuals. GAUDI explicitly models ancestry-differential effects while borrowing information across segments with shared ancestry in admixed genomes. We demonstrate marked advantages of GAUDI over other methods through comprehensive simulation and real data analyses for traits with associated variants exhibiting ancestral-differential effects. Leveraging data from the Women's Health Initiative study, we show that GAUDI improves PRS prediction of white blood cell count and C-reactive protein in African Americans by > 64% compared to alternative methods, and even outperforms PRS-CSx with large European GWAS for some scenarios. We believe GAUDI will be a valuable tool to mitigate disparities in PRS performance in admixed individuals.

Proteomics to predict relapse in patients with myelodysplastic neoplasms undergoing allogeneic hematopoietic cell transplantation.

Disease relapse remains a major barrier to success after allogeneic hematopoietic cell transplantation (allo-HCT) in myelodysplastic neoplasms (MDS). While certain high risk genomic alterations are associated with increased risk of relapse, there is a lack of clinically applicable tools to analyze the downstream cellular events that are associated with relapse. We hypothesized that unique proteomic signatures in MDS patients undergoing allo-HCT could serve as a tool to understand this aspect and predict relapse. Using the Center for International Blood and Marrow Transplant Research (CIBMTR) database, we identified 52 MDS patients who underwent allo-HCT and analyzed their proteomic profile from pretransplant blood samples in a matched case-control design. Twenty-six patients without disease relapse after allo-HCT (controls) were matched with 26 patients who experienced relapse (cases). Proteomics assessment was conducted using the Slow Off-rate Modified Aptamers (SOMAmer) based assay. In gene set enrichment analysis, we noted that expression in the hallmark complement, and hallmark allograft rejection pathways were statistically enriched among patients who had disease relapse post-transplant. In addition, correlation analyses showed that methylation array probes in cis- and transcription regulatory elements of immune pathway genes were modulated and differentially sensitize the immune response. These findings suggest that proteomic analysis could serve as a novel tool for prediction of relapse after allo-HCT in MDS.

Determinants of mosaic chromosomal alteration fitness.

Clonal hematopoiesis (CH) is characterized by the acquisition of a somatic mutation in a hematopoietic stem cell that results in a clonal expansion. These driver mutations can be single nucleotide variants in cancer driver genes or larger structural rearrangements called mosaic chromosomal alterations (mCAs). The factors that influence the variations in mCA fitness and ultimately result in different clonal expansion rates are not well-understood. We used the Passenger-Approximated Clonal Expansion Rate (PACER) method to estimate clonal expansion rate for 6,381 individuals in the NHLBI TOPMed cohort with gain, loss, and copy-neutral loss of heterozygosity mCAs. Our estimates of mCA fitness were correlated (R 2 = 0.49) with an alternative approach that estimated fitness of mCAs in the UK Biobank using a theoretical probability distribution. Individuals with lymphoid-associated mCAs had a significantly higher white blood cell count and faster clonal expansion rate. In a cross-sectional analysis, genome-wide association study of estimates of mCA expansion rate identified TCL1A , NRIP1 , and TERT locus variants as modulators of mCA clonal expansion rate.

Deleterious heteroplasmic mitochondrial mutations are associated with an increased risk of overall and cancer-specific mortality.

Mitochondria carry their own circular genome and disruption of the mitochondrial genome is associated with various aging-related diseases. Unlike the nuclear genome, mitochondrial DNA (mtDNA) can be present at 1000 s to 10,000 s copies in somatic cells and variants may exist in a state of heteroplasmy, where only a fraction of the DNA molecules harbors a particular variant. We quantify mtDNA heteroplasmy in 194,871 participants in the UK Biobank and find that heteroplasmy is associated with a 1.5-fold increased risk of all-cause mortality. Additionally, we functionally characterize mtDNA single nucleotide variants (SNVs) using a constraint-based score, mitochondrial local constraint score sum (MSS) and find it associated with all-cause mortality, and with the prevalence and incidence of cancer and cancer-related mortality, particularly leukemia. These results indicate that mitochondria may have a functional role in certain cancers, and mitochondrial heteroplasmic SNVs may serve as a prognostic marker for cancer, especially for leukemia.

Whole Genome Sequencing Based Analysis of Inflammation Biomarkers in the Trans-Omics for Precision Medicine (TOPMed) Consortium.

Inflammation biomarkers can provide valuable insight into the role of inflammatory processes in many diseases and conditions. Sequencing based analyses of such biomarkers can also serve as an exemplar of the genetic architecture of quantitative traits. To evaluate the biological insight, which can be provided by a multi-ancestry, whole-genome based association study, we performed a comprehensive analysis of 21 inflammation biomarkers from up to 38,465 individuals with whole-genome sequencing from the Trans-Omics for Precision Medicine (TOPMed) program. We identified 22 distinct single-variant associations across 6 traits - E-selectin, intercellular adhesion molecule 1, interleukin-6, lipoprotein-associated phospholipase A2 activity and mass, and P-selectin - that remained significant after conditioning on previously identified associations for these inflammatory biomarkers. We further expanded upon known biomarker associations by pairing the single-variant analysis with a rare variant set-based analysis that further identified 19 significant rare variant set-based associations with 5 traits. These signals were distinct from both significant single variant association signals within TOPMed and genetic signals observed in prior studies, demonstrating the complementary value of performing both single and rare variant analyses when analyzing quantitative traits. We also confirm several previously reported signals from semi-quantitative proteomics platforms. Many of these signals demonstrate the extensive allelic heterogeneity and ancestry-differentiated variant-trait associations common for inflammation biomarkers, a characteristic we hypothesize will be increasingly observed with well-powered, large-scale analyses of complex traits.

A General Approach to Adjusting Genetic Studies for Assortative Mating.

The effects of assortative mating (AM) on estimates from genetic studies has been receiving increasing attention in recent years. We extend existing AM theory to more general models of sorting and conclude that correct theory-based AM adjustments require knowledge of complicated, unknown historical sorting patterns. We propose a simple, general-purpose approach using polygenic indexes (PGIs). Our approach can estimate the fraction of genetic variance and genetic correlation that is driven by AM. Our approach is less effective when applied to Mendelian randomization (MR) studies for two reasons: AM can induce a form of selection bias in MR studies that remains after our adjustment; and, in the MR context, the adjustment is particularly sensitive to PGI estimation error. Using data from the UK Biobank, we find that AM inflates genetic correlation estimates between health traits and education by 14% on average. Our results suggest caution in interpreting genetic correlations or MR estimates for traits subject to AM.

Genetic examination of hematological parameters in SARS-CoV-2 infection and COVID-19.

People hospitalized with COVID-19 often exhibit altered hematological traits associated with disease prognosis (e.g., lower lymphocyte and platelet counts). We investigated whether inter-individual variability in baseline hematological traits influences risk of acute SARS-CoV-2 infection or progression to severe COVID-19. We report inconsistent associations between blood cell traits with incident SARS-CoV-2 infection and severe COVID-19 in UK Biobank and the Vanderbilt University Medical Center Synthetic Derivative (VUMC SD). Since genetically determined blood cell measures better represent cell abundance across the lifecourse, we also assessed the shared genetic architecture of baseline blood cell traits on COVID-19 related outcomes by Mendelian randomization (MR) analyses. We found significant relationships between COVID-19 severity and mean sphered cell volume after adjusting for multiple testing. However, MR results differed significantly across different freezes of COVID-19 summary statistics and genetic correlation between these traits was modest (0.1), decreasing our confidence in these results. We observed overlapping genetic association signals between other hematological and COVID-19 traits at specific loci such as MAPT and TYK2. In conclusion, we did not find convincing evidence of relationships between the genetic architecture of blood cell traits and either SARS-CoV-2 infection or COVID-19 hospitalization, though we do see evidence of shared signals at specific loci.

Differences in the Circulating Proteome in Individuals with versus without Sickle Cell Trait.

BACKGROUND: Sickle cell trait affects approximately 8% of Black individuals in the United States, along with many other individuals with ancestry from malaria-endemic regions worldwide. While traditionally considered a benign condition, recent evidence suggests that sickle cell trait is associated with lower eGFR and higher risk of kidney diseases, including kidney failure. The mechanisms underlying these associations remain poorly understood. We used proteomic profiling to gain insight into the pathobiology of sickle cell trait. METHODS: We measured proteomics ( N =1285 proteins assayed by Olink Explore) using baseline plasma samples from 592 Black participants with sickle cell trait and 1:1 age-matched Black participants without sickle cell trait from the prospective Women's Health Initiative cohort. Age-adjusted linear regression was used to assess the association between protein levels and sickle cell trait. RESULTS: In age-adjusted models, 35 proteins were significantly associated with sickle cell trait after correction for multiple testing. Several of the sickle cell trait-protein associations were replicated in Black participants from two independent cohorts (Atherosclerosis Risk in Communities study and Jackson Heart Study) assayed using an orthogonal aptamer-based proteomic platform (SomaScan). Many of the validated sickle cell trait-associated proteins are known biomarkers of kidney function or injury ( e.g. , hepatitis A virus cellular receptor 1 [HAVCR1]/kidney injury molecule-1 [KIM-1], uromodulin [UMOD], ephrins), related to red cell physiology or hemolysis (erythropoietin [EPO], heme oxygenase 1 [HMOX1], and α -hemoglobin stabilizing protein) and/or inflammation (fractalkine, C-C motif chemokine ligand 2/monocyte chemoattractant protein-1 [MCP-1], and urokinase plasminogen activator surface receptor [PLAUR]). A protein risk score constructed from the top sickle cell trait-associated biomarkers was associated with incident kidney failure among those with sickle cell trait during Women's Health Initiative follow-up (odds ratio, 1.32; 95% confidence interval, 1.10 to 1.58). CONCLUSIONS: We identified and replicated the association of sickle cell trait with a number of plasma proteins related to hemolysis, kidney injury, and inflammation.

The Relationship of Duffy Gene Polymorphism, High Sensitivity C-Reactive Protein, and Long-term Outcomes.

Background: Black adults have higher incidence of all-cause death and worse cardiovascular outcomes when compared to other populations. The Duffy chemokine receptor is not expressed in a large majority of Black adults and the clinical implications of this are unclear. Methods: Here, we investigated the relationship of Duffy receptor status, high-sensitivity C-reactive protein (hs-CRP), and long-term cardiovascular outcomes in Black members of two contemporary, longitudinal cohort studies (the Jackson Heart Study and Multi-Ethnic Study of Atherosclerosis). Data on 4,307 Black participants (2,942 Duffy null and 1,365 Duffy receptor positive, as defined using Single Nucleotide Polymorphism (SNP) rs2814778) were included in this analysis. Results: Duffy null was not independently associated with elevated levels of serum hs-CRP levels once conditioning for known CRP locus alleles in linkage disequilibrium with the Duffy gene. Duffy null status was not found to be independently associated with higher incidence of all-cause mortality or secondary outcomes after adjusting for possible confounders in Black participants. Conclusions: These findings suggest that increased levels of hs-CRP found in Duffy null individuals is due to co-inheritance of CRP alleles known to influence circulating levels hs-CRP and that Duffy null status was not associated with worse adverse outcomes over the follow-up period in this cohort of well-balanced Black participants.

A genome-wide gene-environment interaction study of breast cancer risk for women of European ancestry.

BACKGROUND: Genome-wide studies of gene-environment interactions (G×E) may identify variants associated with disease risk in conjunction with lifestyle/environmental exposures. We conducted a genome-wide G×E analysis of ~ 7.6 million common variants and seven lifestyle/environmental risk factors for breast cancer risk overall and for estrogen receptor positive (ER +) breast cancer. METHODS: Analyses were conducted using 72,285 breast cancer cases and 80,354 controls of European ancestry from the Breast Cancer Association Consortium. Gene-environment interactions were evaluated using standard unconditional logistic regression models and likelihood ratio tests for breast cancer risk overall and for ER + breast cancer. Bayesian False Discovery Probability was employed to assess the noteworthiness of each SNP-risk factor pairs. RESULTS: Assuming a 1 × 10-5 prior probability of a true association for each SNP-risk factor pairs and a Bayesian False Discovery Probability < 15%, we identified two independent SNP-risk factor pairs: rs80018847(9p13)-LINGO2 and adult height in association with overall breast cancer risk (ORint = 0.94, 95% CI 0.92-0.96), and rs4770552(13q12)-SPATA13 and age at menarche for ER + breast cancer risk (ORint = 0.91, 95% CI 0.88-0.94). CONCLUSIONS: Overall, the contribution of G×E interactions to the heritability of breast cancer is very small. At the population level, multiplicative G×E interactions do not make an important contribution to risk prediction in breast cancer.

Whole genome analysis of plasma fibrinogen reveals population-differentiated genetic regulators with putative liver roles.

Genetic studies have identified numerous regions associated with plasma fibrinogen levels in Europeans, yet missing heritability and limited inclusion of non-Europeans necessitates further studies with improved power and sensitivity. Compared with array-based genotyping, whole genome sequencing (WGS) data provides better coverage of the genome and better representation of non-European variants. To better understand the genetic landscape regulating plasma fibrinogen levels, we meta-analyzed WGS data from the NHLBI's Trans-Omics for Precision Medicine (TOPMed) program (n=32,572), with array-based genotype data from the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium (n=131,340) imputed to the TOPMed or Haplotype Reference Consortium panel. We identified 18 loci that have not been identified in prior genetic studies of fibrinogen. Of these, four are driven by common variants of small effect with reported MAF at least 10% higher in African populations. Three ( SERPINA1, ZFP36L2 , and TLR10) signals contain predicted deleterious missense variants. Two loci, SOCS3 and HPN , each harbor two conditionally distinct, non-coding variants. The gene region encoding the protein chain subunits ( FGG;FGB;FGA ), contains 7 distinct signals, including one novel signal driven by rs28577061, a variant common (MAF=0.180) in African reference panels but extremely rare (MAF=0.008) in Europeans. Through phenome-wide association studies in the VA Million Veteran Program, we found associations between fibrinogen polygenic risk scores and thrombotic and inflammatory disease phenotypes, including an association with gout. Our findings demonstrate the utility of WGS to augment genetic discovery in diverse populations and offer new insights for putative mechanisms of fibrinogen regulation. Key Points: Largest and most diverse genetic study of plasma fibrinogen identifies 54 regions (18 novel), housing 69 conditionally distinct variants (20 novel).Sufficient power achieved to identify signal driven by African population variant.Links to (1) liver enzyme, blood cell and lipid genetic signals, (2) liver regulatory elements, and (3) thrombotic and inflammatory disease.

ENIGMA CHEK2gether Project: A Comprehensive Study Identifies Functionally Impaired CHEK2 Germline Missense Variants Associated with Increased Breast Cancer Risk.

PURPOSE: Germline pathogenic variants in CHEK2 confer moderately elevated breast cancer risk (odds ratio, OR ∼ 2.5), qualifying carriers for enhanced breast cancer screening. Besides pathogenic variants, dozens of missense CHEK2 variants of uncertain significance (VUS) have been identified, hampering the clinical utility of germline genetic testing (GGT). EXPERIMENTAL DESIGN: We collected 460 CHEK2 missense VUS identified by the ENIGMA consortium in 15 countries. Their functional characterization was performed using CHEK2-complementation assays quantifying KAP1 phosphorylation and CHK2 autophosphorylation in human RPE1-CHEK2-knockout cells. Concordant results in both functional assays were used to categorize CHEK2 VUS from 12 ENIGMA case-control datasets, including 73,048 female patients with breast cancer and 88,658 ethnicity-matched controls. RESULTS: A total of 430/460 VUS were successfully analyzed, of which 340 (79.1%) were concordant in both functional assays and categorized as functionally impaired (N = 102), functionally intermediate (N = 12), or functionally wild-type (WT)-like (N = 226). We then examined their association with breast cancer risk in the case-control analysis. The OR and 95% CI (confidence intervals) for carriers of functionally impaired, intermediate, and WT-like variants were 2.83 (95% CI, 2.35-3.41), 1.57 (95% CI, 1.41-1.75), and 1.19 (95% CI, 1.08-1.31), respectively. The meta-analysis of population-specific datasets showed similar results. CONCLUSIONS: We determined the functional consequences for the majority of CHEK2 missense VUS found in patients with breast cancer (3,660/4,436; 82.5%). Carriers of functionally impaired missense variants accounted for 0.5% of patients with breast cancer and were associated with a moderate risk similar to that of truncating CHEK2 variants. In contrast, 2.2% of all patients with breast cancer carried functionally wild-type/intermediate missense variants with no clinically relevant breast cancer risk in heterozygous carriers.

On asymptotic distributions of several test statistics for familial relatedness in linear mixed models.

In this study, the asymptotic distributions of the likelihood ratio test (LRT), the restricted likelihood ratio test (RLRT), the F and the sequence kernel association test (SKAT) statistics for testing an additive effect of the expected familial relatedness (FR) in a linear mixed model are examined based on an eigenvalue approach. First, the covariance structure for modeling the FR effect in a LMM is presented. Then, the multiplicity of eigenvalues for the log-likelihood and restricted log-likelihood is established under a replicate family setting and extended to a more general replicate family setting (GRFS) as well. After that, the asymptotic null distributions of LRT, RLRT, F and SKAT statistics under GRFS are derived. The asymptotic null distribution of SKAT for testing genetic rare variants is also constructed. In addition, a simple formula for sample size calculation is provided based on the restricted maximum likelihood estimate of the effect size for the expected FR. Finally, a power comparison of these test statistics on hypothesis test of the expected FR effect is made via simulation. The four test statistics are also applied to a data set from the UK Biobank.

Comparison of multiple imputation and other methods for the analysis of imputed genotypes.

BACKGROUND: Analysis of imputed genotypes is an important and routine component of genome-wide association studies and the increasing size of imputation reference panels has facilitated the ability to impute and test low-frequency variants for associations. In the context of genotype imputation, the true genotype is unknown and genotypes are inferred with uncertainty using statistical models. Here, we present a novel method for integrating imputation uncertainty into statistical association tests using a fully conditional multiple imputation (MI) approach which is implemented using the Substantive Model Compatible Fully Conditional Specification (SMCFCS). We compared the performance of this method to an unconditional MI and two additional approaches that have been shown to demonstrate excellent performance: regression with dosages and a mixture of regression models (MRM). RESULTS: Our simulations considered a range of allele frequencies and imputation qualities based on data from the UK Biobank. We found that the unconditional MI was computationally costly and overly conservative across a wide range of settings. Analyzing data with Dosage, MRM, or MI SMCFCS resulted in greater power, including for low frequency variants, compared to unconditional MI while effectively controlling type I error rates. MRM andl MI SMCFCS are both more computationally intensive then using Dosage. CONCLUSIONS: The unconditional MI approach for association testing is overly conservative and we do not recommend its use in the context of imputed genotypes. Given its performance, speed, and ease of implementation, we recommend using Dosage for imputed genotypes with MAF [Formula: see text] 0.001 and Rsq [Formula: see text] 0.3.

Whole-genome sequencing identifies novel predictors for hematopoietic cell transplant outcomes for patients with myelodysplastic syndrome: a CIBMTR study.

Recurrent mutations in TP53, RAS pathway and JAK2 genes were shown to be highly prognostic of allogeneic hematopoietic cell transplant (alloHCT) outcomes in myelodysplastic syndromes (MDS). However, a significant proportion of MDS patients has no such mutations. Whole-genome sequencing (WGS) empowers the discovery of novel prognostic genetic alterations. We conducted WGS on pre-alloHCT whole-blood samples from 494 MDS patients. To nominate genomic candidates and subgroups that are associated with overall survival, we ran genome-wide association tests via gene-based, sliding window and cluster-based multivariate proportional hazard models. We used a random survival forest (RSF) model with build-in cross-validation to develop a prognostic model from identified genomic candidates and subgroups, patient-, disease- and HCT-related clinical factors. Twelve novel regions and three molecular signatures were identified with significant associations to overall survival. Mutations in two novel genes, CHD1 and DDX11, demonstrated a negative impact on survival in AML/MDS and lymphoid cancer data from the Cancer Genome Atlas (TCGA). From unsupervised clustering of recurrent genomic alterations, genomic subgroup with TP53/del5q is characterized with the significant association to inferior overall survival and replicated by an independent dataset. From supervised clustering of all genomic variants, more molecular signatures related to myeloid malignancies are characterized from supervised clustering, including Fc-receptor FCGRs, catenin complex CDHs and B-cell receptor regulators MTUS2/RFTN1. The RSF model with genomic candidates and subgroups, and clinical variables achieved superior performance compared to models that included only clinical variables.

Prognostic landscape of mitochondrial genome in myelodysplastic syndrome after stem-cell transplantation.

Despite mitochondrial DNA (mtDNA) mutations are common events in cancer, their global frequency and clinical impact have not been comprehensively characterized in patients with myelodysplastic neoplasia (also known as myelodysplastic syndromes, MDS). Here we performed whole-genome sequencing (WGS) on samples obtained before allogenic hematopoietic cell transplantation (allo-HCT) from 494 patients with MDS who were enrolled in the Center for International Blood and Marrow Transplant Research. We evaluated the impact of mtDNA mutations on transplantation outcomes, including overall survival (OS), relapse, relapse-free survival (RFS), and transplant-related mortality (TRM). A random survival forest algorithm was applied to evaluate the prognostic performance of models that include mtDNA mutations alone and combined with MDS- and HCT-related clinical factors. A total of 2666 mtDNA mutations were identified, including 411 potential pathogenic variants. We found that overall, an increased number of mtDNA mutations was associated with inferior transplantation outcomes. Mutations in several frequently mutated mtDNA genes (e.g., MT-CYB and MT-ND5) were identified as independent predictors of OS, RFS, relapse and/or TRM after allo-HCT. Integration of mtDNA mutations into the models based on the Revised International Prognostic Scores (IPSS-R) and clinical factors related to MDS and allo-HCT could capture more prognostic information and significantly improve the prognostic stratification efforts. Our study represents the first WGS effort in MDS receiving allo-HCT and shows that there may be clinical utility of mtDNA variants to predict allo-HCT outcomes in combination with more standard clinical parameters.

Structural variation across 138,134 samples in the TOPMed consortium.

Ever larger Structural Variant (SV) catalogs highlighting the diversity within and between populations help researchers better understand the links between SVs and disease. The identification of SVs from DNA sequence data is non-trivial and requires a balance between comprehensiveness and precision. Here we present a catalog of 355,667 SVs (59.34% novel) across autosomes and the X chromosome (50bp+) from 138,134 individuals in the diverse TOPMed consortium. We describe our methodologies for SV inference resulting in high variant quality and >90% allele concordance compared to long-read de-novo assemblies of well-characterized control samples. We demonstrate utility through significant associations between SVs and important various cardio-metabolic and hematologic traits. We have identified 690 SV hotspots and deserts and those that potentially impact the regulation of medically relevant genes. This catalog characterizes SVs across multiple populations and will serve as a valuable tool to understand the impact of SV on disease development and progression.